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ABSTRACT 

We used the Karl G. Jansky Very Large Array (VLA) to image one primary beam area at 3 GHz with 
8^' FWHM resolution and 1.0 /iJy beam" rms noise near the pointing center. The P{D) distribution 
from the central 10 arcmin of this confusion-limited image constrains the count of discrete sources 
in the 1 < S{fiJy) < 10 range. At this level the brightness-weighted differential count S'^n{S) is 
converging rapidly, as predicted by evolutionary models in which the faintest radio sources are star- 
forming galaxies; and ^ 96% of the background originating in galaxies has been resolved into discrete 
sources. About 63% of the radio background is produced by AGNs, and the remaining 37% comes 
from star-forming galaxies that obey the far-infrared (FIR) / radio correlation and account for most 
of the FIR background at A ~ 160 /im. Our new data confirm that radio sources powered by AGNs 
and star formation evolve at about the same rate, a result consistent with AGN feedback and the 
rough The confusion at centimeter wavelengths is low enough that neither the planned SKA nor its 
pathfinder ASKAP EMU survey should be confusion limited, and the ultimate source detection limit 
imposed by "natural" confusion is < 0.01 /iJy at = 1.4 GHz. If discrete sources dominate the 
bright extragalactic background reported by ARCADE 2 at 3.3 GHz, they cannot be located in or 
near galaxies and most are < 0.03 /iJy at 1.4 GHz. 

Subject headings: cosmology: diffuse radiation — cosmology: observations — galaxies: statistics — 
radio continuum: galaxies 



1. INTRODUCTION 

1.1. Counts of Extragalactic Radio Sources 

Together the number per steradian n{S)dS of discrete 
extragalactic radio sources having flux densities S to 
S -\- dS and the brightness temperature Tb that sources 
contribute to the sky background constrain the nature 
and evolution of all extragalactic radio sources, even 
those too faint to be counted individually. Counts of ra- 
dio sources achie ved instant f ame and notoriety with the 
announcement by Ryle ( 1955 ) and Ryle & Scheuer ( 1955 ) 
that the source count from the 2 (J survey was dramati- 
cally steeper than that expected for a uniformly-filled Eu- 
clidean universe. The remarkable Ryle & Scheuer (1955) 



paper claimed that the majority of "radio stars" were ex- 
tragalactic, that they were so distant and radio-luminous 
as to be mostly beyond the reach of then-existing opti- 
cal telescopes, and that the steep count showed dramatic 
cosmic evolution for individual objects in space density 
or luminosity. Then revolutionary, it remains a succinct 



^ The NRAO is a facility of the National Science Foundation 
operated by Associated Universities, Inc. 



summary of our current understanding of extragalactic 
radio sources. 

The claimed evolution sparked a bitter and public dis- 
pute j_^efutmg_a^Jtd^ cosmol- 
ogy ( jBondi fc Gold|1948l[Hoyle|1948D . Also the 2C radio 
source count was niade before radio astronomers appre- 
ciated the importance of confusion and Eddington bias. 
Mills fc Slee (1957) showed the great majority of the 2C 
sources were just "confusion bumps," the broad 2C beam 
often blending two or more faint back fflound sou r ces to 
appear as a single stronger source. Mills & Slee| (|1957) 



also noted that the count from their higher-resolution 
survey showed "no cosmological effects" with the pos- 
sible exception of mild clustering. By that time, how- 
ever, the Cambridge group had fully understood what 
went wrong, and the aftermath triggered (1) the 3C and 
4C Cambridge surveys which were carefully cognizant of 
co nfusion, ( 2) the pioneering P{D) analysis of confusion 
by Scheuer] ( |1957p which showed how to extract the true 
source count from a confusion-li mited survey, and (3) 
the first aperture-synthesis images dRyle fc Neville|1962p 
which revealed a decreased source-count slope (loosely 
termed "turnover") at flux densities below 1 Jy. 
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By 1966 there was general agreement about the counts 
of sources stronger than S ^ 0.1 Jy at lo w frequen- 
cies. The 4C confusion P{D) (Hewish 1961") and direct 
( Gower||1966 ) counts had mapped their main features: a 
rise steeper than the static Euclidean slope followed by a 
decrease to sub-Euclidean slopes at flux densities below 
1 Jy. Independent sky surveys, such as the Parkes sur- 
vey ( Bolton et al.|p!964| ), were in agreement. To explain 
the source count, several investigators developed mathe- 
matical (non-physical) models for the cosmic evolution of 
extragalactic radio sources. Longair (1966) first showed 
that the relatively rapid change of the count slope could 
be explained by differential cosmic evolution, in which 
the more luminous sources underwent more evolution in 
their comoving space density than their less-luminous 
counterparts did. 

Sensitive surveys made with the WSRT and VLA in 
the 1980s extended the 1.4 GHz source cou nt down to 
mJy levels and below. Condon & Mitchell (fl984) and 



Windhorst et al. (1985) discovered a point oiintlection 
near 1 mJy below which the count slope flattens, suggest- 
ing the emergence of a new population of radio sources. 
Spiral galaxies dominate the local radio luminosity func- 
tion below L ~ 10^^ W Hz~^. Radio sources powered 
by star formation in spiral galaxies can account for this 
new population and for nearly half of the extragalactic 
sky background if they evolve at about the same rate 
as the stronger sources powered by active galactic nu- 
clei (AG Ns) in elliptical galaxies | Condon 1984aL Later 
models (|Wan||1997l iBoyle & Terlevich 1998: mTi 
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. . Boyle & Terlevich 1998. WThnan_et 
evolution similar to that of the cosmic 
star-formation rate shown in the Lilly-Madau diagram 
see [H opkins (2007) for a modern data-set]. Massardi 
et al. l |TQTQp accounted for the inflection with an evolu- 
tionary model using the empirical FIR/radio correaltion, 
FIR counts, and the redshift distributions of star-forming 
galaxies . A recent compilation of source counts from sur- 
veys at frequaencies 1 50 MHz to 15 GHz is presented by 
De Zotti et al.[ (2010); features and cosmic implications 
are discussed there m detail. 

Of particular note in the present context are persistent 
discrepancies among the most sensitive 1.4 GHz direct 
counts t hat far exceed both Poisson and clu stering uncer- 
tainties ("Con don[[2007| lOwen fc Morrison[ r2008). These 
discrepancies may result from ditterent authors making 
different corrections for the effects of partial resolution. 
Resolution corrections can be large and difficult to esti- 
mate for an image whose resolution approaches the me- 
dian angular size (0^) ~ 1" of faint sources. The deep 
synthesis counts were also extended b y confusion P(D) 
distributions in low-resolution ima ges (Mitchell & Con- 
dQii[|1985[ [Windhorst et al.[|1993[ ). We used their ap- 
proach here to count sources a factor of ten fainter in a 
low-resolution (8'' FWHM) 3 GHz image that does not 
need significant corrections for partial resolution. 



1.2. The Source Contribution Sky Brightness 

The radio source count and the source contribution to 
the sky brightness at frequency u are connected by the 
Rayleigh- Jeans approximation 



where /cb ~ 1.38 x 10~^^ J is the Boltzmann con- 
stant and dT\y is the sky brightness temperature added 
by the n{S) sources per steradian with flux densities S 
to S dS. The observed source count now spans eight 
decades of flux density, so it must be plotted with a loga- 
rithmic abscissa. Substituting (i[ln(S')] = dS/S gives the 
brightness temperature contribution per decade of flux 
density 
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The cumulative brightness temperature from all 
sources stronger than any detection limit is 



Arb(> So) 
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SMS)d[\og{S)] . (3) 



If So is low enough, ATb approaches the total Tb of 
the source background, and we can say that the back- 
ground jms_been_^esolved_jnto known sources. For ex- 
ample. 
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(2010) resolved about half of the 
far- infrared (l^lKj background at A = 250, 350, and 
500 /im; and t he Herschel/PEP (PACS Evolutionary 
Probe) survey ([Berta et al.[|2011| resolved about 3/4 of 
the COBE/DIRBE background at A 160 /im. 

Figure [l] shows the flux-density range covered by pub- 
lished 1.4 GHz source counts. We plotted the weighted 
count log[6'^n(5')] instead of the traditional Euclidean- 
hted count log[S^^'^n{S)] as a function of log(5') be- 
cause (1) S'^n{S) is proportional to the source contri- 
bution per decade of flux density to the sky temper- 
ature (Eq. [2|, (2) the radio universe is neither static 
nor Euclidean, and (3) the plotted slopes are minimized 
for easy visual recognition of broad features, unlike the 
steeply sloped plot of log[n(5')] for example. On a plot 
of log[5'^n(6')], the source count must ultimately fall off 
at both ends to avoid Olbers' paradox. 

The fllled points at log\S (Jy)] > — 3 are from the many 
surveys referenced in[Condon[ (|1984a[) , and the fllled data 
points at lower flux densities are from the Mitchell 



Condon ( [l985) confusion-limited VLA survey. The ir 
regular box encloses the range of 1.4 GHz counts consis- 
tent with the probablility distribution P{D) of confusion 
amplitudes D (in brightness units /iJy per beam solid 
angle) in the low-resolution (17V5 FWHM) Mi tchell fc 
Condon (1985) image. The straight line inside the box 
indicates the best power-law flt to the P{D) data. The 
open data points and their power-law flt (upper straight 
line) indicate the direct count of individual sources from 
the most sensitive high-resolution 1.4 GHz survey ever 
made with the original VLA (Owen fc Morrison 2008). 
These two faint-source counts disagree by several times 
the Poisson counting errors, and nearly all other pub- 
lished faint-source counts are scattered over the wide gap 
b etween them ( Condon 2007^, as illustrated by Figure 11 



m 



Owen fc Morrison ( 2008 ) . The sources contributing to 



Sn{S)dS = 
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the data plotted in Figure [T] show no sign of resolving the 
extragalactic background at 1.4 GHz because the power- 
law flts to S'^n{S) are still rising near the flux density 
limit ^0 ~ 10 /iJy. 

The solid curve in Figure [l] shows the [Condon[ ( [1984b[ ) 
evolutionary model for the total source count at 1.4 GHz. 
In this model, the stronger radio sources are powered 



Resolving the Radio Source Background 



3 



primarily by AGNs of elliptical galaxies (dashed curve) 
and most of the fainter sources are in star-forming spi- 
ral galaxies (dotted curve). The model predicts that (1) 
S'^n{S) should converge below 5 ~ 10 fiJy and (2) most 
of the model's Tb ^ 100 mK background should be re- 
solved into sources stronger tha n 5* ~ 1 /iJy. 

However, Fixsen et al. (2011) reported that the recent 
ARCADE 2 (Absolute Radionieter for Cosmology, As- 
trophysics, and Diffuse Emission) balloon measurement 
of absolute sky brightness leaves a remarkably high ex- 
tragalactic Tb = 54 ± 6mK at u = 3.3 GHz after the 
Galactic foreground and the T = 2.731 ± 0.004 K cosmic 
microwave background (CMB) were subtracted. Com- 
bining their 3.3 G Hz result with low- f requency data from 
the literature led Fixsen et al.| ( [2QlT ) to fit the "excess" 
extragalactic background with a power-law spectrum 



Tb = (24.1 ±2.1 K) X 



310 MHz 



-2.599±0.036 



(4) 



between 22 MHz and 10 GHz. Equation [4] gives Tb ^ 
480 mK at 1.4 GHz, almost five times the Tb ^ 100 mK 
predicted for the known populations of extragalactic 
sources. 

With the goals of (1) determining an accurate source 
count down to S ^ 1 /iJy, (2) resolving the radio sky 
background contributed by galaxies, and (3) constrain- 
ing possible source populations that could produce the 
high ARCADE 2 background at 3.3 GHz, we used the 
Karl G. Jansky Very Large Array (VLA) to make a very 
sensitive (rms noise dn ~ l/iJybeam~^) low-resolution 
(8^^ FWHM) confusion-limited sky image covering one 
primary beam area at S band (2-4 GHz). We pointed 
at the center of the "crowded" Owen & Morrison ( 2008 ) 
field to see if we could confirm the reported high source 
count there. Our relatively low angular resolution en- 
sures that galaxy-size sources at cosmological distances 
are unresolved: 8^' spans at least 40 kpc throughout the 
redshift range OA < z < 7 containing most faint radio 
sources. That low resolution also guarantees that the 
rms confusion is at least comparable with the rms noise, 
a necessary condition for using confusion to constrain 
the sky density of sources as faint as ~ 1 /iJy, which 
is more than a factor of five below our detection limit for 
individual sources. 

1.3. Outline 

Section |2] of this paper describes our observations and 
the production of a confusion-limited sky image. Sec- 
tion [3] presents the 3.02 GHz P{D) distribution from the 
central region of this image. The 1.4 and 3.02 GHz source 
counts consistent with our P{D) distribution are derived 
in Section |4| and the contributions of known extragalac- 
tic source populations to the sky brightness are discussed 
in Sectionp] In Section [6 we use our narrow P{D) distri- 
bution to snow that thenigh ARCADE 2 background is 
too smooth to be produced by, or even spatially associ- 
ated with, galaxies brighter than ttiab = +29. Section [7| 
summarizes our results. 

2. OBSERVATIONS AND IMAGING 

2.1. Target Field 



over the 2-4 GHz S-band frequency range using an 
average of 21 antennas in the C configuration (maxi- 
mum baseline length ^ 3 km) of the VLA. Th is field in 
the Locknian Hole was originally selected by |Owen fc 
Morrison] (|2008) because it was covered by the Spitzer 
Wid e- area InfraR ed Extragalactic (SWIRE) legacy sur- 
vey ( Lonsdale et al.|2003 ), has exceptionally low infrared 
cirrus, and contains no strong radio sources. We fa- 
vored it because it is far enough north to minimize S- 
band radio-frequency interference (RFI) from powerful 
geosynchronous broadcast satellites just south of the ce- 
lestial equator and has the highest reported faint-source 
count (Owen fc Morrison^^2008J . 

2.2. Observations 

Our 57 hours of observing time was divided among six 
observing nights between 2012 February 21 and March 
18, and 50 hours of this time was spent integrating on 
the target field. The nearby unresolved phase calibra- 
tor J1035+564 was monitored for 30 seconds every 30 
minutes. The flux density and bandpass calibrators 3C 
147 and 3C 286 were each observed once per night. An- 
tenna pointing accuracy was maintained with corrections 
determined from J1035+564 every hour. Accurate point- 
ing is necessary even at S band because pointing errors 
can easily limit the dynamic range of an image filled with 
sources (Condon 2009). 

The data- averaging period was 1 second of time, and 
the 2-4 GHz frequency range was divided into 16 corre- 
lator subbands, each with 64 spectral channels of width 
2 MHz. The spectral resolution was broadened and 
the averaging intervals were increased somewhat prior 
to imaging, but both bandwidth and time smearing were 
always kept small enough that our sky coverage is limited 
only by the primary-beam attenuation of the individual 
25 m antennas. Their attenuation pattern is very nearly 
that of a unifor mly illuminated circular aperture (Con- 
don et al.|[T998l ): 



2Ji(3.233p/6>) 
(3.233p/6>) 



(5) 



where Ji is the Bessel function of the first kind and order, 
p is the angular offset from the pointing center, and 



l9p = (43f3 ± 0;4) 



GHz 



(6) 



is the measured full width between half-maximum points 
(FWHM) of the S-band primary beam at frequency u. 
The primary FWHM ranges from 21 J 7 at u = 2 GHz to 
10.^8 at u = 4GHz. 

2.3. Editing and Calibration 
We used the Obit package (Cotton 2008 Q to edit and 



We reobserved the "crowded" Owen fc Morrison ( 2008 ) 



field centered on J2000 a = 10^ 46"^ 00^ ^ = +59° 01^ 00^' 



calibrate the {u^ v) data. The six observing sessions were 
calibrated and edited separately. Two of the 16 sub- 
bands contained satellite RFI strong enough to cause se- 
rious Gibbs ringing in the raw data, a consequence of 
the finite number of correlator lags. Hanning smooth- 
ing (combining adjacent spectral channels with weights 
1/4, 1/2, and 1/4) suppressed this ringing and broadened 

^ http://www.cv.nrao.edu/~bcotton/Obit.html 
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the spectral resolution to about 4 MHz. The new VLA 
correlator was designed to prevent "bleeding" of ringing 
and other nonlinearities from one subband to another. 
Spectral channels still containing strong interfering sig- 
nals were flagged and removed from the data, as were a 
few edge channels in each subband. Prior to calibration, 
the remaining data samples containing strong interfer- 
ing signals that were impulsive in time or frequency were 
identified by their large deviations from running medi- 
ans in time and frequency, and they were fiagged. Cali- 
brator observations were also subjected to a test of the 
RMS/mean amplitude ratios in each data stream, and 
noisy data segments with high RMS / mean amplitude ra- 
tios were deleted. 

Subsequent calibration and editing consisted of the fol- 
lowing steps: 

1. Instrumental group delay offsets were deter- 
mined from observations of 3C 147, 3C 286, and 
J1035+564 and applied to all of the data. 

2. Residual variations of gain and phase with fre- 
quency were corrected by bandpass calibration 
based on 3C 286. 

3. Amplitude calibration was based on the VLA stan- 
dard spectrum of 3C 286 (see Table 1, column 3) 
bootstrapped to determine the spectrum of the 
astrometric calibrator J1035+564. The more fre- 
quent observations of J1035-I-564 were then used 
to calibrate the amplitudes and phases the target 
(ix, v) data. Data from some antennas, time in- 
tervals, and frequency ranges were still degraded 
by interference that had evaded earlier editing and 
corrupted a small fraction of our amplitude and 
phase solutions. These corrupted solutions were 
detected by a comparison with all solutions, and 
data implying deviant complex gains were fiagged. 

4. The calibrated data were then subjected to further 
editing in which data with excessive Stokes I or V 
amplitudes were deleted, as well as another pass at 
removing narrowband interference. 

After this calibration and editing, the initial calibra- 
tion was reset and the whole process was repeated us- 
ing only the data that had survived the editing pro- 
cess. Finally, because the target field contains no sources 
brighter than 5mJy beam~^, we were able to fiag the 
small amount of the data having amplitudes significantly 
above the noise. About 53% of the (li, v) data survived 
all of the editing steps. The calibrated and edited (ix, v) 
data from all six observing sessions were then combined 
for imaging. 

2.4. Imaging 

Observations spanning the frequency range between 
the low frequency limit ui and the high frequency limit 
z^h have center frequency Uc = (z^i + z^h)/2, bandwidth 
Au = (z/h — ^i), and fractional bandwidth Au/uc The 
fractional bandwidth covered by our 2-4 GHz (ix, v) data 
is exceptionally large: Avjv^ = 2/3. Using such data 
to make an image suitable for measuring confusion en- 
counters three interesting problems — the field of view, 
the point-spread function (PSF), and the fiux densities 



of most sources can vary significantly with frequency. To 
deal with these problems, we separated the (u^v) data 
into 16 subbands each having a small fractional band- 
width Avjv^ <C 1. The [u^v) data were tapered heav- 
ily in the higher-frequency subbands, and each s ubband 
was imaged with an independent "robustness" (Briggs 
1995[ ) to force nearly identical PSFs in all subbands. We 
weighted and recombined the narrowband images to pro- 
duce a sensitive wideband image characterized by an "ef- 
fective frequency" (z^), where {v) at any position in a 
wideband image is defined as the frequency at which the 
fiux density of a point source with a typical spectral in- 
dex (a) = dXviS I dXviV = —0.7 equals its fiux density 
in the wideband sky image. The effective frequency de- 
clines with angular distance p from the pointing center 
because the primary beamwidth is inversely proportional 
to frequency. 

Table 1 lists the center frequencies of the 16 sub- 
bands and the rms noise values dn of the 16 subband 
images. The fractional bandwidths of these subbands 
range from 3% to 6%, so each subband image is a nar- 
rowband image. We used the Obit task MFImage to form 
separate dirty and residual images in each subband. We 
assigned each subband image a weight inversely propor- 
tional to its rms noise, generated a combined wideband 
image from their weighted average, and used this sensi- 
tive combined image to locate CLEAN components for 
a joint deconvolution. The CLEAN operation used fiux 
densities from the individual subband images but at lo- 
cations selected from the combined image. At the end of 
each major CLEAN cycle, the CLEAN components with 
fiux densities from the individual subband images were 
used to create residual subband {u^ v) data and then new 
residual images. 

The effects of sky curvature across the field of view 
were eliminated by faceted "fiy's eye" imaging in which 
each facet was small enough to make the curvature ef- 
fects negligible. Each facet image was projected onto a 
common coordinate grid on one tangent plane to form 
a single continuous central image 30 arcmin in radius. 
This allowed all facets to be CLEANed in parallel. Out- 
side the central image, separate ou tlier facets were added 
at the positions of those NVSS (Condon et al.||1998) 



sources whose uncleaned sidelobes might otherwise at- 
feet our target field. Before imaging, the data were av- 
eraged over baseline-dependent time intervals subject to 
the constraints that (1) the averaging should not cause 
time smearing within our central image and (2) the av- 
eraging time should never exceed the 20 second phase 
self-calibration interval. 

In order to obtain nearly identical PSFs in all sub- 
bands, we tapered the (ix, v) data in each subband differ- 
ently and assigned individual "robust" weighting factors 
adjusted to ensure that each synthesized "dirty" beam 
was nearly circular with major and minor axes between 
7" and S" . After CLEANing, each residual image was 
smoothed by convolution with its own elliptical Gaussian 
tailored to yield a circular and nearly Gaussian PSF with 
a precisely 8^0 FWHM. Having the same dirty beam and 
PSF in each subband image is critical for an accurate con- 
fusion analysis because only those sources much stronger 
than the rms confusion were actually CLEANed. Finally, 
all CLEAN components were restored with an 8^' FWHM 
circular Gaussian beam. 
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Two iterations of phase-only self calibration were used 
to remove residual atmospheric phase fluctuations. The 
final CLEAN cycle was stopped when the peak in the 
combined residual image reached 10/iJy beam~^, a level 
well above the rms noise and confusion in the combined 
image. The total CLEANed fiux density in the combined 
image is 21 mJy. 

Unwanted fiuctuations indistinguishable from confu- 
sion are produced by any dirty-beam sidelobes remain- 
ing in our combined image. Fortunately, they are small 
because the combination of long observing tracks and 
bandwidth synthesis over our wide fractional bandwidth 
ensures excellent (ii, -plane coverage and keeps the 
dirty-beam sidelobe levels well below 1% of the peak, 
as shown in Figure [2] Consequently the highest dirty- 
beam sidelobes from sources in the residual image are 
< 0.01 X 10/iJy beam~^ ^ 0.1 /iJybeam"^, so their con- 
tribution to the image variance is more than two orders of 
magnitude below the ^ (1 A^Jy beam~^)^ contributions 
of noise and confusion. 

We used the AIPS task IMEAN to calculate the rms 
noise values an of the CLEANed subband images in 
several large areas that are well outside the main lobe 
of the primary beam and contain no visible (5'p > 
6/iJybeam~^) sources; the are listed in Table 1. 
Next we assigned to each subband image a weight in- 
versely proportional to its noise variance <j^, generated 
a wideband image from the weighted average of the sub- 
band images, and measured the noise distribution in four 
large regions well outside the primary main beam and 
free of visible sources. Figure [3] shows a logarithmic his- 
togram of the noise amplitudes in one such region cover- 
ing an area = NVti, containing N ^ 4000 synthesized 
beam solid angles Vt\^ = 7r^g/(41n2) ^ 0.020 arcmin^ for 
= The excellent parabolic fit to the logarithmic 
histogram indicates a precisely Gaussian noise amplitude 
distribution with an rms ~ 1.02 /iJy beam~^. 

What is the rms statistical uncertainty Adn in our esti- 
mate of (Jn? The image PSF is a Gaussian, and the effec- 
tive noise area for a Gaussian PS F is t hat of the Gaus- 
sian PSF squared (Condon 19971 [Cond o n et al . 1998). 
Squaring a Gaussian PSF of widtn 6 yields a narrower 
Gaussian of width 2~^/^6> and solid angle r2b/2. Conse- 
quently there are actually 2N statistically independent 
noise samples in ^7 = NQ]j beam solid angles, and the 
rms fractional uncertainty in (Jn is 



AcTn 



1-) 

2N 



1/2 



(7) 



not the commonly beliefed (1/A^)^/^. See the Appendix 
for a detailed derivation of equation [71 

The final result from the four areas covering a to- 
tal N = 11570 beam solid angles is cr^ = 1.012 ± 
0.007 /iJy beam" ^ (statistical error only). In theory the 
rms noise is uniform across the image prior to correction 
for primary-beam attenuation, so we used this value to 
estimate the noise in confusion-limited regions near the 
pointing center. The total intensity-proportional error 
arising from uncertainties in the fiux-density calibration 
and primary beamwidth is not more than 3% inside the 
primary beam half-power circle. 



2.5. The SNR- Optimized Wideband Sky Image 

Our final wideband image of the sky was made with 
weights designed to correct for primary-beam attenu- 
ation and simultaneously maximize the signal-to-noise 
ratio (SNR) for sources having spectral indices near 
(a) = —0.7, the mean spectral index of faint source s 
found at frequencies around 3 GHz (Condon 1984b). 
This differs from the traditional weighting designed to 
minimize noise, which maximizes the SNR in a narrow- 
band image, but in a wideband image only if (a) ~ 0. 
The brightness bi{p) of each pixel in each of the z = 1, 16 
subband images was assigned a weight 



CrnA(p, i/c) 



(8) 



where the rms noise and center frequency Uc of each 
subband image is listed in Table 1. Each pixel in the 
weighted wideband sky image was generated from the 
ratio 

16 / 16 

b{p) = Y.Mp)Wi{p)] / J2w^{p) (9) 

i=l / i=l 

Even at the pointing center, weighting to optimize the 
SNR for (a) = —0.7 increases the sky image noise 
slightly from = 1.012 ± 0.007 /iJy beam~^ to = 
1.080 ± 0.007 /iJy beam~^. Away from the pointing cen- 
ter, the frequency-dependent primary-beam attenuation 
correction causes the rms noise on the sky image to grow 
with the radial offset p at the rate indicated by the 
dashed curve in Figure ffl Weighting also affects how 
the effective frequency (i^/at each point in the wideband 
sky image decreases with the offset p from the pointing 
center. The dashed curve in Figure pi indicates how (u) 
in our sky image declines monotonically from 3.06 GHz 
at the pointing center to 2.96 GHz at p = 5arcmin. 

Figure |6] is a profile plot of the central portion of our 
SNR-optimized wideband sky image. The intensity scale 
can be inferred from the highest peaks, which are trun- 
cated at Sp = 100 /iJy beam~^. This image is confusion 
limited in the sense that the rms fluctuations are every- 
where larger than the noise levels plotted in Figure |4] 

3. THE P{D) DISTRIBUTION 

Historically, confusion was measured from the proba- 
bility distribution P{D) of pen deflections D on a chart- 
recorder plot of fringe amplitudes from the si ngle bas eline 
of a two-element radio interferometer (Scheuer|1957). On 



a single-dish or aperture-synthesis array image, the cor- 
responding "deflection" D at any pixel is the intensity in 
units of flux density per beam solid angle. The observed 
D at any point is the sum of the contribution from the 
noise- free source confusion and the contribution from im- 
age noise. The source confusion and noise contributions 
are independent of each other, so the the observed P{D) 
distribution is the convolution of the source confusion 
and noise distributions, and the variance of the ob- 
served P{D) distribution is the sum of variances of the 
noise- free source confusion (a^) and the noise (a^) dis- 
tributions: 

al=a^,+al. (10) 
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The goal is to extract the source confusion distribution 
and its width 

a. = {al-alf/' (11) 

from the observed deflections and the measured noise. 
If (Jc < CTn, then dcFc/dcFo ^ {cFo/cf^) ^ {cf^/cf^) > 1 
and small errors in the measured <Jo (or a^) cause large 
errors in the extracted cr^. This is the reason that only 
fairly low-resolution (ctc > (Jn) images can be used to 
measure confusion. Also, to the extent that other sources 
of error (e.g., uncleaned sidelobes) are present but not 
accounted for, equation [TT] will tend to overestimate the 
source confusion. 

The noise in our weighted sky image is low at the cen- 
ter but increases with radial distance p from the pointing 
center, and the noise eventually overwhelms the confu- 
sion at large p. On the other hand, increasing the radius 
p of the circular region inside which the P{D) distribu- 
tion is measured increases the number N of beam solid 
angles sampled and thus acts to decrease the statistical 
uncertainty Actq in the width of the observed P{D) dis- 
tribution. For each thin ring of radius p covering N beam 
solid angles, we use equation [lO] to calculate 



(Aao)^ 



dcFp 
dcFr. 



Act, 



where 



Then 



AcTo 



and 



dcFp 
dan 



d(Jn 



AcTr 

_ ^ 

Act, 



(12) 



(13) 



(14) 



where (Aac/ac) A/""^/^ and (Aan/an) 
(Eq. [t]). In the limit <Jn ^ ctc, 

2 



(2Ar)- 



-1/2 



Ao-o 



(15) 



To minimize the uncertainty ctq inside the circle of radius 
p, we counted each pixel a"^ times and truncated the 
count at p ^ 5 arcmin, where AcTo/cTo in the enclosed 
circle has a broad minimum. Within the circle of radius 
p = 5 arcmin, our weighted P{D) distribution has rms 
noise (Jn = 1.255 /iJy beam"^ (Fig. [f]) and its effective 
frequency is (u) ~ 3.02 GHz (Fig. pi). 

The weighted P{D) distribution within 5 arcmin of the 
pointing center is shown by the data points with their 
Poisson rms error bars in Figure ^ The thin curve in- 
dicates the Gaussian noise distribution with rms width 
CTn = 1.255 /iJy beam" ^. Interferometers have no DC 
response to smooth emission, so the mean deflection 
(D) ^ 3/iJy beam~^ in our image approximately equals 
the total CLEANed flux density divided by the number 
of beam solid angles in the image. Consequently, the 
image zero-point D = should be treated as a free pa- 
rameter when comparing the observed P{D) distribution 
with analytic models of the true sky P{D) distribution. 

The noise-free confusion P{D) distribution can be cal- 
culated analytically for the scale-free case of a power-law 
differential source count ( Condon 1974 J 



where k is the count normalization and 1 < 7 < 3 is the 
differential count slope. The shape of the noise-free P{D) 
distribution depends only on 7, and the width of the 
noise-free P{D) distribution obeys the scaling relation 



Piikn^y^^^-^^D] = (me)"^/^^"^^P(L)) , 



(17) 



where Qq is the effective beam solid angle. If G is the 
normalized (G = 1 at the peak) PSF, then 



I 



and 



7- 



(18) 



(19) 



for a circular Gaussian PSF. Equation [Tsl also implies 
that sidelobes of the PSF have a bigger effect as 7 de- 
clines, so images suitable for confusion studies at sub-/iJy 
levels, where 7 — 1 ~ 0.5, need to have very low sidelobe 
levels. The rms confusion contributed by deflections 
fai nter than some limiting signal-to-noise ratio q = Sq/cTc 
is ( |Condon||1974| ) 



^3-7 \ 1/(7-1) 

C7 = I l^^j {kn,)'/^^-'^ 



(20) 



If the number of sources per steradian stronger than S is 
N{>So)= / n{S)dS (21) 

then there are 



/3=[N{> So)n^]-' (22) 

beam solid angles per source stronger than 5*0. For a 
power-law source count, 

P=i^ ' (23) 



■7 



Taking g = 5 for reliable source detection ( [Murdoch et| 
al. 1973) and 7^2, equation 20 limits the maximum 
density of reliably detectable anoindividually countable 
sources to about one source per 25 beam areas if 7 ^ 2. 

These equations apply to unresolved sources. The re- 
sponse to an extended source is the convolution of the 
source brightness distribution with the point-source re- 
sponse, and areas add under convolution. The integrated 
flux density is conserved, so the peak flux density falls. 
In the limit of a slightly resolved source covering a solid 
angle Vt^ <C l^e^ the peak flux density will be multiplied 
by the factor 

/ ^ < 1 (24) 

and the solid angle of the source response on the image 
will be divided by /. Thus equation [l^ can be used to 
show that slightly extended sources multiply the width 
of the P{D) distribution by 



(25) 



n{S) = kS- 



(16) 



y.(7-2)/(7-l) ^ 

The thick curve in Figure [7| is the convolution of the 
<jn = 1.255 /iJy beam~^ noise Gaussian with the noise- 
free P{D) distribution for the power-law source count 
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n{S) = 90005-^-^ Jy"' sr-^ of point sources and an 8'^ 
FWHM Gaussia n PSF. Even if faint sourc es are as large 



as (6>s) ^ l'/5 dOwen fc Morriso^ |2QQ8'), equation [25 



indicates that the width of the i^(Lij3istribution woula 
increase by less than 2%. Thus the source count based 
on our low-resolution P{D) data should be much less 
sensitive to the angular-size distribution of faint sources 
than direct source counts from images that must have 
synthesized beam solid angles about 25 times smaller to 
reach the same flux-density limit. 

4. SOURCE COUNTS AT 1.40 AND 3.02 GHZ 

To compare our 3.02 GHz P{D) distribution with 
source counts at 1.4 GHz, we assume that the average 
spectral index between 1.4 and 3.02 GHz is (a) = —0.7, 
the value that best fits multifrequency counts of stronger 
sources (Condon 1984b). The power-law count kiS~^ at 
z/i = 1.4 GHz can be converted to the count k2S~^ at 
V2 = 3.02 GHz via the relation 



(«>(7-l) 



(26) 



If our estimate of (a) is in error by 0.05 and 7 = 1.7, 
the frequency conversion will cause an ^ 3% error in the 
calculated value of /c2. 
To span the range of faint-source counts reported at 



1.4 GHz, w e converted the relatively low [Mitchell fc Con^ 
(|1985|) co unt n{S) = 



don 



. 57 5-^ - ^ Jy~ ^ sr-^ and the 

relatively nigh |Owen fc Morrison (2008) count n{S) = 
QS~'^-^ Jy~^ sr~^ at ui = 1.4 GHz to counts at U2 = 
3.02 GHz and extrapolated them to lower flux densi- 
ties. The noise-free 3.02 GHz P{D) distributions for 
these counts observed with a ^ = 8^' FWHM Gaussian 
beam were convolved with the = 1.255 /iJy beam~^ 
Gaussian noise distribution for coniparison with the ob- 
served P{D) distribution (Figure [8|). The continuous 
curve is slightly lower and broader than the P(D) data, 
indicating that the extrapolated [Mitchell fc Condon] 
( p85l ) count is somewhat high, as expected fro m the 
rapi d ly conv erging weighted count S'^n{S) of the Con- 



don 



(1984b) model. The dashed curve is much lower 



and br oader because the extrapolated [Owen fc Morrison 



(2008) count is much higher. 

Next we compared our 3.02 GHz P{D) distribution 
with the predictions of two 1.4 GHz ev olution models: 

(1) The very simple ^ondonl (|1984bl) model is based on 
old local 1.4 GHz luminosity functions derived fro m op- 
tically selec ted samples of ellipt ical galaxies (,Auriemma 
et al.[1977D and spiral galaxies ( |Hummel[[T98(jl ) . I'he ra- 
dio sources in elliptical and spiral galaxies are primarily 
associated with AGNs and star formation, respectively, 
so these luminosity functions are similar to more mod- 
ern luminosity functions in which the FIR/radio corre- 
lation and/or optical line spectra are used to distinguish 
radio sources primarily powered by AGN from those pri 



marily powered by recently fo rmed stars (Condon et al. 
|2002^ 'Mauch fc Sadler"2007'). There were only sparse 
redshift data for identifications of strong radio sources, 
and a pre-WMAP cosmology was assumed. The decision 
to evolve both the AGN and star-forming source popu- 
lations identically to match the source count at several 
frequencies was based on applying Occam's razor to lim- 
ited source-count data below S ^ 1 mJy, and not any 



physical theory. 

(2) The fa r more sophisticated and modern Wilman 
et al.[ ( |2008[ ) semi-empirical simulation of the sky was 
made as part of the Square Kilometer Array Design 
Study (SKADS). It uses new and better local luminosity 
functions, redshift data, and a concordance cosmological 
model incorporating dark matter and dark energy. 

Remarkably, both models predict nearly identical 
source counts all the way down to 5* lOnJy, and 
both agree well with our observed P{D) distribution. 
This good agreement suggests that there is little "wig- 
gle room" for models describing the counts of faint 
sources once the local radio luminosity function has been 
specified; the details of the evolution and the underly- 
ing cosmology see m to have littl e effect on the result- 
ing source count (|Condon[[T989|. The fact that our 
P{D) count matches so well the S' ^n(S) con vergence be- 
low S ^ 10/iJy pred icted by the Condon ( |1984b) and 



Wilman et al. (2008) models provides independent sup- 
port tor the claim that most faint radio sources are in 
star-forming galaxies. Such radio sources are usually 
somewhat smaller than their optical host galaxies, so we 
expect that the median angular diameter of /iJy sources 
is Os < I". 

In the flux-density range ^ ^ S < 10 /iJy, the 
1.4 GHz model counts can be approximated b y the 
pow er laws n{S) = 1.2 x lO^^"^-^ Jy"^ sr-^ ( |Con- 
d^ [l984b|) and n(S) = 3.45 x 10"^ 6' 



sr 

-1-6 Jy 



sr 



( | Wilman et ar|[2008D . At 3.02 GHz these counts become 
n{S) = 9.17 X 10^-1-^ Jy"^ sr-^ and n{S) = 2.5 x 
lO'^^'-^-^ Jy"^ sr- \ respect i vely, if (a) = -0.7. Fig- 



ure [9| compares the 



Condon[ ( [T984b| (dashed curve) and 



Wilman et al. (2008) (continuous curve) model P{D) dis- 
tri but ions, after conversion to 3.02 GHz with (a) = —0.7 
and convolution with = 1.255 /iJy beam~^ noise, with 
the 3.02 GHz P{D) data. 

Although the power-law approximations to both mod- 
els appear to agree fairly well with the P{D) data in 
Figure |9j the normalized xl values from fits between 



D 



-3.9 /iJy beam and D = +9.7/iJybeam 



= 66 degrees of freedom) are too high: 



for the power-law approximation to the Condon ( 1984b ) 
m odel and ^1 for the power-law approximation to the 



100 



Wilman et al.[ (2008) model. Both approximations sig- 
niticantly overestimate the numbers of sources stronger 
than about 5 jaJy. The P{D) data are simply not con- 
sistent with any power-law approximation to the source 
count over the whole range 

1 ^ ^(^Jy beam"^) < 10. 
The best-fit slope at the high end of this range is 7 ^ 1.7 
while the best-fit slope at the low end is 7 ~ 1.5. 
Better modeling of the source counts will require nu- 
merical simulations to generate P{D) distributions for 
non-power-law source counts. The range of 3.02 GHz 
power-law counts reasonably consistent with the P{D) 
data is spanned by the three fits alr eady mentioned : 



n{S) = 9.17 X lO'^S-^-^ Jy"^ sr "^ ([Condon 
n{S) = 2.5 X lO^S-^-^ Jy"^ sr-^ ( [Wilman et al. 
and n{S) = dOOOS'^-^ Jy"^ sr-^ {Fig.^. 



1984b 



2008 



4.1. The Noise-Free RMS Confusion 
The corresponding noise-free P{D) distributions for a 



= 8'' FWHM Gaussian PSF are plotted in Figure 10 
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Note that the zero-point deflection D = in these cal- 
culated P{D) distributions is the absolute zero of sky 
brightness contributed by discrete sources. For values of 
7 < 1.7, all sources fainter than ^ 1 /iJy contribute 
very little to the sky background. The rms confusion dc 
is a poor statistic for describing such skewed P{D) dis- 
tributions. In a Gaussian distribution, about 2/3 of the 
points lie within ila of the mean, so we use a* defined 
as half the range of D containing 2/3 of the points as a 
more stable measure of the width of the confusion dis- 
tribution. In a i9 = 8'' FWHM PSF, this width of the 
noiseless confusion distribution is a* ^ 1.2 /iJy beam~^ 
at u = 3.02 GHz. For nearby values of and the 
scaling equation [20] with 7 ^ 1.6 implies 



cr* ~ 1.2 /iJy beam ^ 



3.02 GHz 



-0.7 



10/3 



• (2T) 



For example, the proposed EMU (Extragalactic Map 
of the Universe) survey will cover most of the sky at 
u = 1.3 GHz with = W resolution ( |Norris et al. 12011 ) 
and rms noise <Jn = 10 /iJy beam"^.^ Equation |27| in- 
dicates that EMU will not be confusion limited, with 
crj ^ 5 /iJy beam~^. 

4.2. The Natural Confusion Limit 

In addition to the instrumental confusion limit caused 
by finite angular resolution, there is a "natural" con- 
fusion limit caused by the finite angular size of faint 
sources. If the median solid angle covered by a faint 
source is (r^s), then equations 



22 



and 



23 imply the num- 



ber of source solid angles per oetectaHe source {q > 5) 
must exceed 



(3 = [N{> s)n,] = 



3-7 



25 



(28) 



even if the instrumental resol ution is very high. Their 
high source count prompted lOwen fc Morrison (2008|) 
to speculate that the natural contusion limit could be 
as high as 1 /iJy, but our lower P{D) count suggests 
that the natural confusion limit for detecting individ- 
ual (fls) ~ (I'O^ sources is under 10 nJy. This is well 
below the hoped-for = 45 nJy sensitivity limit of the 
fuh 12000 m^ K-^ Square Kilometre Array (SKA) af- 
ter 1 megasecond of integration, so the SKA continuum 
sensitivity should not be limited by either natural or in- 
strumental confusion, but only by instrumental noise and 
dynamic range (Condon 2009). 

4.3. The 14 GHz Source Count 

The full range of 1.4 GHz source count data is plot- 
ted with the traditional S^^'^n{S) "static Euclidean" nor- 
malization in Figure [TT] for direct comparison with most 
published results. The range of 1.4 GHz source counts 
consistent with our P{D) distribution is conservatively 
within the box covering 2 < /S < 20 /iJy. Again, ob- 
taining more precise constraints on the counts consistent 
with our P{D) distribution will require numerical simu- 
lations of non-power-law counts. 

The upper abscissa of Figure pT] shows the source- frame 
1.4 GHz spectral luminosity L(W Hz~^) of a normal- 
spectrum {a ^ —0.7) source whose fiux density is di- 
rectly below on the lower abscissa if the source is at the 



median re dshift (z) ~ 0.8 of extragalactic sources (Con- 
don||1989" ). These typical luminosities and flux densities 



at z ^ ir8 can be compared with those of nearby radio 
sources such a s the Large M agellanic Cloud (LMC) at 
log(L) ^ 20.3 ( [Hughes et aL| |2007) and \ok (S) ^ - 7.1, 
our own Galaxy at log(L) « 20 (Berkhuijsenll 984) and 
\og(S) ^ -6.0, M82 at log(L) ^ 22.1 (Keherma nn et al. 
|1969| ) and \o^(S) ^ -5.3, and Arp 220 at log(L) « 23.4 
( jCondon et al.|2002[ ) and log(5') ^ -4.0. Most of the ex- 



tragalactic radio sources fainter than our P{D) limit were 
no more luminous at z ~ 0.8 than our Galaxy is today. 
For the first time, th e count data are deep enough to con- 
firm the assumption ( Condon 1984a|b ) that radio sources 
powered by AGNs and radio sources powered by star for- 
mation evolve at about the same rate. That assumption 
had no known physical justification in 1984, but it is 
consistent with the correlation between black hole and 
stellar bulge masses ( ,Magorrian et al. 1998) discovered 
later. Note that this correlation appears to be weaker for 
the lower-mass (^ 10^ M^) black holes fo und in late-type 
star- forming galaxies ( Greene et al.|201Q ). The processes 
that produced the correlation tor niassive elliptical galx- 
ies do not appear to have had enough time to fully es- 
tablish the correlation for the star-forming galaxies that 
host most faint radio sources. 

Our new source count is about a factor of four lower 
than the jOwen fc Morrison) ( |2008 ) count near 6* = 15 /iJy, 
even though our tields overlap on the sky. What might 
cause this discrepancy and, for that matter, the sur- 
prisingly large scatter among all published faint-source 
counts jCo ndon 2007)? Are most of the faint [Owen 



Morrisoui (^2008j) sources spurious , or did we miss a 
large traction of real sources that |Owen fc Morrison 



( |2008[ ) counted? The Owe n fc Morrison[ (120081 ) image 
has 6 = 1 '.'6 angular resolu tion, much~mgn er than our 



resolution and the Mitchell fc Condon (^1985 



6 = 17 V5 resolution. How often have we blended into 



one "source" what Owen fc Morrison ( 2008 ) resolved into 
two or more sources? I^'igure|l2|shows our 3 GHz sky im- 
age as a gray scale with white crosses on the positions 
of cataloged Owen fc Morris'onl (|2008) sources stronger 
than 15 /iJy at 1.4 GHz! The correspondence is actu- 
ally quite good. Almost all crosses have 3 GHz coun- 
terparts, and most of the brighter 3 GHz sources have 
1.4 GHz counterparts, so spurious or missing sources can- 
not explain a factor of four discrep ancy. Likewise , there 
are only a few cases of two Owen fc Morri soni (|2 008[ ) 
sources blending into one source on our low-resolution 
image. We suspect that most of the count difference 
is caused by count corrections made for partial resolu- 
tion of extended sources in the high-resolution 1.4 GHz 
beam. Survey catalogs are complete to a fixed bright- 
ness (/iJy beam~^) cutoff, so extended sources with lower 
brightnesses but higher integrated fiux densities (/iJy) 
will be missed. The corrections needed to convert source 
brightnesses in /iJy beam~^ to source fiux densities in 
/iJy and to account for missing sources become quite 
large near the brightness cutoff as the angular resolu- 
tion approac hes the median angu l ar siz e of faint sources. 
For example, Owen fc Morrison (2008) found a median 
angular size {(p) = 1'.'5 in their 20-30 /iJy bin, which 
is about the same as their FWHM resolution = IV6. 
The rapid rise in counts near cutoffs at 6 , 7, 8, 9, and 10 
times the rms noise shown in Figure 9 of Owen fc Morri- 
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son ( 2008 ) may be a sign that these corrections were too 
large. 

5. THE CONTRIBUTIONS OF KNOWN SOURCES TO THE 
SKY BRIGHTNESS 

The 1.4 GHz source counts in Figure [TT] were replotted 



as log[6' n(6')] versus log(6') in Figure 13 to emphasize 



the source contribution to sky brightness. The right or- 
dinate of the top panel in Figure [13] shows the differential 
contribution dT]j / d[\og{S)] to the 1.4 GHz sky brightness 
temperature per decade of flux density. The constraint 
provided by the new P{D) distribution clearly shows that 
S'^n{S) is falling off at low flux densities, at about the 
rate predicted by the Con don] ( |1984b[ ) and jWilman et al.| 
( [2008, ) evolutionary models. 

We estimated the 1.4 GHz and 3.02 GHz backgrounds 
from the known source populations powered by AGNs 
and star- forming galaxies by inserting the 1.4 GHz source 
counts indicated by the solid curve in the top panel of 
Figure 13 into equation [3] and assuming an effective spec- 
tral index (a) = —0.7. The results are shown in the 
middle and bottom panels, respectively. The total back- 
ground from the model sources converges to Tb « 100 mK 
at 1.4 GHz and Tb ^ 13 mK at 3.02 GHz. AGNs account 
for about 63% and star-forming galaxies 37% of the ex- 
tragalactic source background. Approximately 96% of 
this total has been resolved into sources stronger than 
So ^2 /iJy at 1.4 GHz or S'o 1 jaJy at 3.02 GHz. More 
specifically, nearly 100% of the AGN background and 
about 89% of the star-forming galaxy contribution to the 
background has been resolved. 

Star-forming galaxies appear to obey the local 
FIR/radio correlat ion even at moderate redshifts z ^ 2 
( Ivison et al.||2010" ), while the radio sources powered by 
AGINs have much lower FIR/radio flux ratios. Local 
{z <^ 1) star- forming galaxies have a mean FIR/radio 
flux ratio (q) = log [5(80 /im)/6'( 1.4 GHz)] ^ 2.3 (Con- 



|don[ 

duce 



T992]). 
ITTb 



Our claim that star-forming galaxies pro- 
^ 37 mK background at = 1.4 GHz 
(equivalently, - 2.2 x 10"^^ W m"^ Hz"^ sr"^) 



X iU W m Mz sr 
can be checked by comparison with the measu r ed ex- 
tragalactic 
and [Lagache et al. 



FIR bac k ground . 

|T999) 



Fixsen et al. 

analyzed COBF 

13.7 ± 4 nW m"^ sr"^ at 



(1998) 
FIRXS 



data to find uljy 

160 /im, or ^ 7.3 x 10"^^ W m"^ Hz"' sr-^ at 
u ^ 1870 GHz. Thus the background FIR/radio ra- 
tio is (7b = log[/^(160/im)//^(1.4 GHz)] ^ 2.5. If 
the star-forming background galaxies have typical red- 
shifts (z) ~ 0.8 and radio spectral indices a ~ —0.7, 
then these backgrounds imply that their source-frame 
FIR/radio ratio is log[^(90 /im)/^(2.5 GHz)] ^ 2.5, 
which corresponds to a source- frame FIR/radio flux ra- 
tio log[6'(90/im)/6'(1.4 GHz)] ^ 2.3. This agrees with 
the local FIR/radio correlation and supports the idea 
that the same population of galaxies accounts for the 
star-forming radio and FIR backgrounds. 

The Herschel PACS Evolutionary Probe (PEP) re- 
solved 74 ± 5 % of the A = 160 /im background into in- 
divi dually detec t ed g alaxies stronger than about 3 mJy, 
and |Berta et al. (2011) used a P{D) analysis to raise the 
resolved traction to ^ 89% above ~ 0.3 mJy. The sensi- 
tivity limits of the VLA {0 = 8'' FWHM resolution) and 
Herschel PEP A = 160 /im {0 ^ 11" FWHM resolution) 
P{D) distributions to galaxies obeying the FIR/radio 



correlation are about equal, and the same fraction of 
the star-forming galaxy background has been resolved at 
radio and FIR wavelengths. Converting the direct and 
P{D) counts at A = 160 /im to 1.4 GHz using the re- 
lation log[6'(160/im)/S'(1.4 GHz)] = 2.5 allows them to 
be compared directly with the Con don, (,1984 b) 1.4 GHz 
count model for star-forming galaxies (dotted curve in 
Fig. 14); the agreement is well within the statistical er- 
rors. 

6. THE ARCADE 2 EXTRAGALACTIC BACKGROUND AND 
CONSTRAINTS ON NEW POPULATIONS OF RADIO 
SOURCES 

If the extragalactic source background is really Tb ~ 
100 mK at 1.4 GHz and 13 mK at 3.02 GHz, then it has 
been resolved into the known source populations powered 
by AGNs and star formation. However, all of the faint- 
source counts are based on interferometric images which 
are insensitive to a smooth background, the 2.7 K CMB 
for example. We cannot rule out a sufficiently smooth 
distribution of faint radio sources that contributes signif- 
icantly to the extragalactic background brightness, but 
we can set strong lower limits to the numbers of sources 
needed to make the background smooth enough to agree 
with our narrow 3.02 GHz P{D) distribution. 

The ARCADE 2 measurement of the extragala ctic sky 
temperatur e in the 3-90 GHz frequency range (| Fixsen| 
et al. 2011|) recently reported a surprising extragalactic 
excess over the CMB. Its brightness spectrum is 



(24.1 ±2.1) 



0.31 GHz 



-2.599±0.036 



(29) 



between 22 MHz to 10 GHz. This is T5 = 480±50mK at 
u = lA GHz and Tb = 65±8 mK at = 3.02 GHz, much 
higher than the lOOilOmK and 13±1.3mK attributable 
to known extragalactic radio sources at 1.4 GHz and 
3.02 GHz, respectively. 

One possible explanation is the existence of a new pop- 
ulation of faint extragalactic sources that contributes 
ATb ^ 480 - 100 = 380 ± 50 mK at 1.4 GHz and 
ATb ^ 65 - 13 = 52 ± 8mK at 3.02 GHz to the sky 
background. A truly diffuse synchrotron background 
produced by relativistic electrons in interg alactic space 
far fro m individual galaxies was ruled out by |Singal et al. 
(2010) on the grounds that inverse- Compton scattering 
of those electrons radiating in weak {B < 1 fiG) inter- 
galactic magnetic fields would exceed the ob served X-ray 
and 7-ray backgrounds. Singal et al.| (2010') also pointed 
out that the diffuse emission from clusters of galaxies 
has a spectrum that is too steep (a < — 1) to explain the 
ARCADE 2 excess, as do es emission from shocks in the 

the source count of any new pop- 
sources contributing ATb to 




the sky brightness. Figure 15 plots log[5'^n(6')] as a func 



tion of log(5') for the known populations of radio sources 
powered by AGN and by star- forming galaxies; they ap- 
pear as two broad peaks centered near log(6'pk) ~ — 1 
and log(6'pk) ~ — 5, respectively. The regions near these 
peaks contribute most to the sky brightness, and each 
peak is well represented by the parabolic Taylor-series 



approximation log[6' n(5')] 
S^n{S) ^ Aexp\-4\n{2) 



-6[log(5)-log(Vr or 



[log(5)-log(V)]^ 



}■ 



(30) 
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where (j) is the logarithmic FWHM of the Gaussian and 
iSpk is the flux density of the peak. The peak amplitude 
A (the peak brightness in units of Jy sr~^) required for a 
new population characterized by a Gaussian of FWHM (j) 
to add ATb to the background can be found by inserting 
equation 30 into equation [3] and integrating over all flux 
densities (6o = 0); it is 



2 r 



Act): 



ln(10)c2 



ln(2) 



nl/2 



(31) 



A and (j) are individually free, but their product A(j) is 
fixed by ATb. 

Figure [15] shows three examples of source counts, each 
of which yield ATb ^ 380 mK at = 1.4 GHz. The dot- 
ted, dashed, and continuous curves correspond to loga- 
rithmic FWHMs ^ = 0.2, 1.0, and 2.0 in equation [3Q1 and 
their peak amplitudes calculated from equation |31| are 
log[A(Jysr"^)] ^ 4.67, 3.97, and 3.67, respectively. The 
value of ATb is independent of 5'pk, so these parabolic 
curves are free to shift horizontally when constrained by 
sky brightness alone. However, a large value of de- 
creases the number of sources per square arcmin needed 
to produce the excess background and increases the Pois- 
son fluctuations in sky brightness, so the width of our ob- 
served P{D) distribution sets upper limits to the values 
of 5'pk associated with these three examples. If the back- 
ground is too bright and smooth, the number of faint 
radio sources required will exceed the total number of 
galaxies. 

The narrowest possible peak (^ <C 1) needs the small- 
est number N of sources per beam solid angle to pro- 
duce a given Tb. In that unrealistic limit, all of the new 
sources have nearly the same very low flux density iS'pk 
and their noiseless P{D) distribution is nearly Gaussian. 
To set a very conservative upper limit to the rms width of 
this Gaussian, we found the widest Gaussian consistent 
with our observed P{D) distribution. It is shown in Fig- 
and its rms is a ~ 0.70 /iJy beam~^ ^ 1.47 mK 
3.02 GHz. To avoid excessive Poisson noise, 
the minimum number N of sources per beam solid angle 
lib ^ 0.020 arcmin^ must be 



ure 
at 



16 



ATb 



52 mK 
1.47 mK 



1.3 X 10^ 



(32) 



implying > 6 x 10^ radio sources per arcmin^, and the 
maximum source flux density consistent with the back- 
ground brightness is 5'pk < 34nJy at 1.4 GHz if a = 
—0.7. This is slightly below the hoped-for ba^ = 45nJy 
sensitivity limit of the fuh 12000 m^ K-^ SKA after 1 
megasecond of integration, so even if such faint sources 
exist, it will be a long time before they can be detected 
individually. 

Broader flux distributions can only raise the minimum 
source density and lower the peak flux density. If ^ = 1 
((/) = 2), then there must be more than 1.6 x 10^ (3 x 
10^) sources per arcmin^ with peaks at 5'pk < 22nJy 
(< 5nJy). The dotted, dashed, and continuous curves 



in Figure 15 are as far right as allowed by the P{D) 



smoothness constraint. They could be shifted to the left, 
but that would only increase the implied sky density of 
radio sources. Our new constraint on the smoothness of 
the background forces the "bump" in the weighted source 



count at /iJy levels propo sed by Seiffert et al. (2011 ) and 
Vernstrom & Wall (2011) to much lower flux levels. 

There are only 10"^ galaxies brighter than ttiab = +29 
in the 11 arcmin^ Hubble Ultra Deep Field (HUDF) 
(Beckwith et al.| |2006| , or ^ 10^ galaxies per arcmin^, 
about two orders magnitude lower than the minimum 
sky density of faint radio sources needed to produce the 
smooth ARCADE 2 excess background. The HUDF can 
detect galaxies 2 mag fainter than T* out to redshifts 
z ~ 6, so it appears that the smooth ARCADE 2 back- 
ground cannot be produced by galaxies or by objects lo- 
cated in or near individual galaxies (e.g., multiple radio 
supernovae or gamma-ray bursts in galaxies, radio-quiet 
quasars, radio halos smaller than 8 arcsec in diameter) or 
by star-forming galaxies that are more radio-loud than 
expec ted from the FIR/radio correlation (S ingal et al.| 
2010). Might faint sources produced by WIM F annihi- 
lations or d ecays in dark-matter (DM) halos (Fornengo 
et al. 2011) be sufliciently numerous to produce a suf- 
ticiently sniooth background? Annihilation in DM ha- 
los around galaxies would violate our smoothness con- 
straint, leaving only the possibility of preferential DM 
annihilation in high-redshift mini- ha los. If our 3.02 GHz 
P{D) smoothness constraint and the Fixsen et al. (2011) 



3.3 GHz excess brightness are both correct, then a very 
numerous {N > 10^^ over the whole sky) and unexpected 
population of radio sources not associated with galaxies 
has been discovered. 

7. SUMMARY 

We used 21 antennas in the Karl G. Jansky VLA 
to observe a single fleld at S band (2-4 GHz) with a 
FWHM resolution = 8" and reached an rms noise 
(Jn ~ 1 /iJy beam~^ near the image center after 50 hours 
of integration time. The image is confusion limited 
with an "rms" confusion level a* ~ 1.2 /iJy beam~^ at 
u = 3.02 GHz. 

The 3.02 GHz differential source count was derived 
from the confusion P{D) distribution. For comparison 
with published source counts, we converted it to 1.4 GHz 
via the effective spectral index {a) ^ —0.7. The power- 
law approximation n{S) ex to the 1.4 GHz source 
count has a slope approaching 7 ^ 1.5 near 1 /iJy, signif- 
icantly lower than the slopes of published counts above 
10 /iJy. If the faintest radio sources have median angu- 
lar size (^s) ^ f as expected for emission coextensive 
with star-forming regions in distant galaxies, the natu- 
ral confusion limit for source detection is not more than 
5(7* ^ 0.01 /iJy at 1.4 GHz, and the continuum sensitiv- 
ity of the planned SKA will not be limited by natural 
confusion. The observed count is well flt by evolutionary 
models in which the local radio luminosity functions of 
all sources associated with both AGNs and star forma- 
tion evolve at the same rate. This is broadly consistent 
with the correlation of black hole and stellar bulge masses 
in massive elliptical galaxies, although this correlation is 
not yet well established for the lower-mass black h oles 
and bulges in late- type galaxies ( Greene et al.|[2QTQ ). 

The brightness-weighted count S"^n{S) is clearly con- 
verging below 10 /iJy. Our image has resolved about 
96% of the radio background produced by all galax- 
ies (Tb ^ 100 mK at 1.4 GHz and Tb ^ 13 mK at 
3 GHz). Nearly 100% of the ^ 63 mK AGN-powered 
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background at 1.4 GHz has been resolved. The remain- 
ing ^ 37 mK comes from star-forming galaxies that obey 
the FIR/radio correlation and account for most of the 
extragalactic background at A = 160 /im. We resolved 
about 89% of the star-forming galaxy contribution. 

The ARCADE 2 balloon experiment indicated a non- 
thermal excess brightness over the Galaxy, the CMB, and 
that expected from known populations of radio sources in 



galaxies. At 3.02 GHz this excess brightness temperature 
is 52 ± 8 mK. Our narrow 3.02 GHz P{D) distribution 
implies that the excess background must be very smooth. 
Any new discrete-source population able to produce such 
a bright and smooth background is far too numerous to 
be associated with galaxies brighter than ttiab = +29. 
Facilities: VLA 



APPENDIX 

NOISE DISTRIBUTIONS IN SYNTHESIS IMAGES 

In the noise-limited outer regions of our image, the pixels have a Gaussian brightness distribution with rms dn ~ 
1 /iJy beam~^. What is the rms fractional uncertainty (Aan/cr^) of the rms noise averaged over a large solid angle 
Q ^ The answer is complicated because the pixel values are not independent. Just as the limited (u^v) plane 
coverage of data going into our image acts as a convolving filter to smooth our image of the sky brightness distribution, 
it acts as a filter to smooth the image noise distribution. The smoothing function for both sky brightness and noise 
is the synthesized beam. Our synthesized beam is closely approximated by a circular Gaussian (Fig. [2]) with FWHM 
Os = 8 arcsec, so its normalized (peak gain G = 1) gain profile can be written 



G((9, 0) = exp 
The corresponding synthesized beam solid angle is 



-4 In 2 



(Al) 



/z]2 
Gdn = ^ « 1.13^2 (A2) 

It is easy to calculate (A<Jn/crn) for a square top-hat beam with FWHM W: 

G{x, y) = l, - W/2 <x< +1^/2, - W/2 <y< (A3) 

and G{x^y) = elsewhere, so the beam solid angle is Q\) = W^. A large area Q can be tiled hy N = ^ 1 

nonoverlapping beam solid angles. Averaging their N independent noise values divides the rms noise by A/"^/^, so 



in this case. This sort of argument seems to underly the common but incorrect belief that equation A4 is true for all 
beamshapes. 

The problem is that noise variance <j^, rather than noise rms dn, adds under convolution. The amount of noise 
smoothing is proportional to the "noise beam solid angle" 



/ 



= / G'dn (A5) 



rather than to the beam solid angle given by equation |A2| Only in the case of a top-hat beam does = ^b- The 
square of a Gaussian is also a Gaussian, but its FHWM is d^j^}!^ so l^n = for a Gaussian beam. Thus 

. /P, X 1/2 / 1 \ 1/2 

Aa. (^.\ ( ^\ (A6) 



Vl ) \2N 

To convince skeptics that this result correctly describes our data, we selected 31 indepenent source-free regions from 
our image, uncorrected for primary attenuation so dn ~ 1 /iJy beam~^ in all regions. Each region covers 100 x 100 = 10^ 



pixels ~ 1^252 on a side, so Q\) ~ 46.3 pixels and there are N ^ 216 ^ 1 beam solid angles per region. If equation A4 
were correct for our Gaussian synthesized beam, we would expect 



CTn 

Equation | A6| predicts the lower value 

Act, 



0.068 . (A7) 



0.048 . (A8) 



In the 31 regions the observed mean and its rms scatter are 

A(Jn 



= 0.040 ± 0.007 . (A9) 
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This is consistent with equation A6 



but is four standard deviations below equation A4 



A similar argument applies to calculations of survey speed. How much time is needed to make a survey covering 
a large solid angle ft = N^h using a single Gaussian beam? If staring time r is needed to reach the desired survey 
sensitivity at t he center of one bea m, the time needed to make the whole survey is t ^ 2Nr, twice the frequently 
quoted t ^ Nr ( [Condon et ar]|1998[ ). 

Finally, the effective beam area for confusion (Eq. 18) becomes Qq = f G^dQ = when 7^3. In this limit, 
the confusion fluctuations are dominated by the very faintest sources and the P{D) distribution becomes a Gaussian, 
just like a noise distribution. 
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Figure 1. Published 1.4 GHz source counts. The brightness- weighted source count S'^n(S) is proportio nal to the contribution per decade 
of flux density to the sky background temperature Tb. The fllled points at \og[S (Jy)] > —3 are from |Condori| ( 1 984a) and fMitchell &] 
[Condon ([1^985). The polygon encloses the range of 1.4 GHz counts consistent with confusion (Mit chell Sz Uondonp985 ), and the straight 
line mside the box is the best power-law flt to the confusion data. The open data points and their power-law nt (upper straight line) 
indicate the |Owen Sz Morrisonj ( |200 8 ) source count. The solid curve is the Condon, ( |1984bj ) model count composed of sources powered 
primarily by AGNs (dashed curve) and by star formation (dotted curve). Left ordinate: log of the 1.4 GHz source count S'^n(S) (Jy sr"-*^). 
Right ordinate: log of the source contribution dT^^ / d[\og(S)] (mK) to the 1.4 GHz background per decade of flux density. 
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Figure 2. Dirty-beam profile. The dirty beam in our wideband image, shown truncated at 20% of its peak to demonstrate that the 
highest sidelobe level is under 1%. 
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Figure 3. Noise distribution. This noise distribution from one source- free region near the edge of our wideband image is a nearly perfect 
Gaussian with rms <j = 1.02)LtJy beam~^; it appears as a parabola on this logarithmic plot. Abscissa: Log of the number of pixels per bin 
of width 0.2)LtJy beam~^. Ordinate: Peak flux density (micro Jy /beam). 
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Figure 4. Noise levels in the final image. The rms noise <Tn in the ring of radius p in the sky image (dashed curve) and in the weighted 
P{D) distribution of points inside the circle of radius p. Abscissa: Offset from the pointing center (arcmin). Ordinate: Root-mean-square 
noise (/LtJy beam~^). 
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Figure 5. Effective frequencies in the final image. The effective frequency (u) is the frequency at which the flux density of a source with 
spectral index a = —0.7 equals the local flux density in the ring of radius p in the sky image (dashed curve) or the average flux density 
of the weighted P{D) distribution of points inside the circle of radius p. Abscissa: Offset from the pointing center (arcmin). Ordinate: 
Effective frequency (GHz). 
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Figure 6. Confusion profile. This profile plot shows the 3 GHz confusion amplitude in an 8 arcsec FWHM beam, truncated at 100 
microJy/beam. 
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Figure 7. Observed P(D) distribution. The observed probability distribution P(D) (data points with error bars) of "deflections" or 
peak flux densities at 3.02 GHz within a circle of radius p = 5 arcmin. The PSF is a circular Gaussian with FHWM 6 = 8". The 
noise has a Gaussian amplitude distribution with rms <7 n — 1.255 juJy / beam (thin curve). The thick fltted curve is the convolution of the 
noise distribution with the noise-free P{D) distribution ( |Condon|l974| of the power-law source distribution n{S) = 9000S'~-'^-^ Jy""*^ sr"-*^. 
Abscissa: Deflection or peak flux density (/iJy/beam) at 3.02 GHz. Ordinate: Probability density (/iJy beam"-*^). 
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Figure 8. Expected P(D) distributions from extrapolations of published counts. This figure compares the observed P(D) distribution 
(data points with error bars) at 3.02 GHz with power-law extrapolations of deep source counts measured at 1.4 GHz, assuming a mean 
spectr al index (a) = — 0-7. T he continuous curve indicates the expected 3.02 GHz P(D) distribution from the power-la w extrapolation of 



the Mitchell & Condon ( 1985| 1.4 GHz source count n{S) = 57S~'^-'^ ^y~^ sr~-^ , and the dashed curve corresponds to the |Owen Morrison 
( |2008| 1.4 GHz source count n{S) = 6S~'^-^ ^y~^ sr~-^. Abscissa: Defiection or peak fiux density (/iJy/beam) at 3.02 GHz. Ordinate 
Probability density (/iJy beam"-*^). 
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Figure 9. Expected P(D) distributions from source-count models. This figure compares the observed P(D) (data points with error bars) 
at 3.02 GHz with power-law fits to model source counts at 1.4 GHz, assuming a mean spectral index (a) = —0.7. The co ntinuous curve 
indicates the expected 3.02 GHz P(D) distribution for the 1.4 GHz source c ount n(S ) = 3.45 x 10^ S*"-*^-^ predicted by Wilman et all 
(2008fc and the dashed curve corresponds to n{S) = 1.2 x lO^S* predicted by | C ondon] ( fl 9 84b J . Abscissa: Defiection or peak fiux density 
(/iJy/beam) at 3.02 GHz. Ordinate: Probability density (/iJy beam"-*^). 
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Figure 10. Noise-free P{D) distributions. The range of noise-free 3.02 GHz P{D) distributions in a FWHM = 8" Gaussian beam calcu- 
lated for three powe r-law approximations to the 3.02 GHz source cou nt that are consistent with our data: n{S) = 250005'" ■'^•^ Jy""*^ sr"-*^ 
( jWilman et al.|2008t , n{S) = OlTOOS'-i-^ Jy"^ sr-^ ( |Condon|l984bt , and n{S) = QOOOS'-^ '^ Jy-^ sr-^ (this paper). Abscissa: Deflection 
or peak flux density (/iJy/beam) at 3.02 GHz. Ordinate: Probabihty density (/xJy beam~^). 
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Figure 11. Euclidean-normalized source counts. The plotted source ource count s at 1.40 GHz are based on individual sources (data 
points), the 1.4 GHz P{D) distribution (blue line and surrounding error box) from 'Mitchell &: Gondon (1985), and our 3.02 GHz P{D) 
distribution (black line and surrounding error box) convert ed to 1-40 GH z via a mean spectral index a = —0.7. The red and green lines 
are power-law approximations to the 1.40 GHz models of G ondon| ( |1984b| > andjWilman et al. ( 2008 ) , respectively. The dashed curve is the 
contribution of AGNs and the dotted curve is the contribution ot star-torming galaxies, trom the G ondon| ( |1984bj > model. Lower abscissa: 
log flux density (Jy) at 1.40 GHz. Upper abscissa: log spectral luminosity (W Hz"-*^) at 1.4 GHz in the source frame for sources at the 
typical redshift (z) ~ 0.8. Ordinate: Source count normalized by the source count n(S) oc S~^^'^ in a static Euclidean universe (Jy^/^ sr"-*^). 
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Figure 12. Gray-scale sky image at 3 GHz. Sources in the [Owen Morrison| ( [2008| > catalog are indicated by crosses on our gray-scale 
image. The wedge indicates 3 GHz brightnesses between —10 and +20 /xJy beam"^ 
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Figure 13. Brightness weighted source counts and model sky background temperatures at 1.4 and 3.02 GHz. The top panel shows the 
weighted 1.40 GHz source count S'^n(S), the weighting that reflects the contribution of sources to the sky background temperature Tb per 
decade of flux density S. The plotted ource counts at 1.40 GH z are based on individual sources (data points), the 1.4 GHz P(D) distribution 
(blue line and surrounding error box) from Mitchell & Condon (1985), and our 3.02 GHz P(D) distribution (black line and surrounding 
error box) converted to 1.40 GHz via a mean spectral mdex a = — U.7. The red and green lines are power-law approximations to the 1.40 
GHz models of Condon (1984b) and Wilman et al. (2008), respectively. The dashed curve is the contribution of AGNs and the dotted curve 
is the contribution ot star-torming galaxies, trom the Uondon (1984b) model. Abscissa: log flux density (Jy). Ordinates, top panel: Log of 
the 1.40 GHz source count S^n{S) (Jy sr"-*^) on the left, weighted to be proportional to the source contribution dT]^ / d[\og{S)] (mK) to the 
sky background per decade of flux density on the right. Ordinate, middle panel: Log of the 1.40 GHz sky background temperature (mK) 
contributed by all sources stronger than S at 1.40 GHz. Ordinate, bottom panel: 3.02 GHz sky background temperature (mK) contributed 
by all sources stronger than S at 3.02 GHz. 
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Figure 1 4. FIR source count converted to 1.4 GHz. The data points and P{D) error box show the A = 160 /xm FIR counts of 'Berta et] 
|al.|(|2011| ) converted to 1.4 GHz via the FIR/radio relation log[S'(160/im)/S'(1.4 GHz)] = 2.5. The curves are predictions from the (JondoiJ 
(|l^84b|> 1.4 GHz count model for sta r-forming galaxies (dots), AGNs (dashes), and their sum (continuous). The agreement with the actuaf 
1.4 (JjJz counts shown in Figure [13] is within the errors. Abscissa: log flux density (Jy). Ordinates: Log of the 1.40 GHz source count 
S^n(S) (Jy sr"-*^) on the left, weighted to be proportional to the source contribution dT\^ / d[\og(S)] (mK) to the sky background per decade 
of flux density on the right. 
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Figure 15. Source counts at 1.4 GHz consistent with the ARCADE 2 background and our P(D) distribution. Source count s at 1.40 GHz 
based on individual sources (data points), the 1.4 GHz P{D) distribution (blue line and surrounding error box) from Mitche ll Co ndon] 
|1985^, and our 3.02 GHz P(D) distribution (black line and surrounding error box) converted to 1.40 GHz via a mean spectral mdex 
a = —0.7. The red and green lines are power-law approximations to the 1.40 GHz models of Condon (1984b) and Wilman et al. (20081, 
respectively. The dashed curve is the contribution of AGNs and the dotted curve is the contribution ot star-forming galaxies, trom tJie 
Ipondon (1984b) model. The three parabolas at nanoJy levels correspond to models for a hypothetical new population of radio sources 
able to produce the ARCADE 2 excess b ack ground temperature. Their widths are (j) = 0.2 (dotted), 1.0 (dashed), and 2.0 (continuous 
curve), where phi is defined by equation |30| Our upper limit on the rms background fluctuation sets upper limits on the flux densities 
Spk cit the peaks of these parabolas and nence lower limits on the sky density of the new population. Even in the limit of a very narrow 
flux-density distribution, the number of sources per square arcmin is N > 6 x 10^ arcmin"^, which greatly exceeds the N ~ 10^ arcmin"^ 
galaxies brighter than m = +29 in the Hubble Ultra Deep Field (HUDF). Abscissa: log 1.4 GHz flux density (Jy). Ordinates: log of the 
1.40 GHz source count S^n(S) (Jy sr"-*^) on the left, weighted to be proportional to the source contribution dT\^ / d[\og{S)] (mK) to the sky 
background per decade of flux density on the right. 
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Figure 16. The largest possible contribution of numerous faint sources to the observed P(D) distribution. A sharp "bump" in the number 
of extremely faint (5 <^ 1 yLtJy) faint sources would contribute a Gaussian to the observed 3.02 GHz P(D) distribution. The rms of that 
Gaussian cannot exceed a ^ 0.70 /LtJy beam" ^ without excessively broadening our observed P{D) distribution, no matter how low the source 
count near S ~ lA^Jy. Abscissa: Deflection or peak flux density (/iJy/beam) at 3.02 GHz. Ordinate: Probability density (/LtJy beam~^. 
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Table 1 

Subband Images 



Subband Frequency S(3C 286) C7n 



number 


(CjrJdzj 


(Jy) 


()LtJy beam 


01 


2.0500 


12.260 


9.220 


02 


2.1780 


11.877 


13.190 


03 


2.3060 


11.522 


18.170 


04 


2.4340 


11.193 


4.476 


05 


2.5620 


10.885 


4.244 


06 


2.6900 


10.597 


4.314 


07 


2.8180 


10.327 


4.039 


08 


2.9460 


10.073 


3.719 


09 


3.0500 


9.877 


4.477 


10 


3.1780 


9.648 


3.145 


11 


3.3060 


9.432 


3.067 


12 


3.4340 


9.226 


2.965 


13 


3.5620 


9.031 


2.882 


14 


3.6900 


8.845 


3.468 


15 


3.8180 


8.668 


4.447 


16 


3.9460 


8.580 


4.402 



